Predictions of the residue cross-sections for the elements Z=113 and Z=114 
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An extremely good reproduction of experimental excitation function of the In reactions producing 
Z=110,Z=111 and Z=112 is obtained by the two-step model and the statistical decay code KEWPIE. 
Thus, an extension of the recipe permits us to predict reliable values of the residue cross-sections of 
the elements Z=113 and Z=114 which will be a useful guide for planning of experiments. 
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PACS numbers: 25.70.Jj, 24.60.-k, 25.60.Pj, 27.90.+b 



INTRODUCTION 

The production and the identification of new elements 
have been made since the discovery of theperiodic ta- 
ble of the elements by Mendelejeff in 1869 Nowdays 
we are still looking for the heaviest elements, so called 
super- heavy elements(SHE). The properties and the con- 
ditions of the existence of the SHE have been foreseen 
since the establishment of the nuclear shell model 
The studies of shell structure predict an island of sta- 
bility for those elements around the double closed shell 
next to ^^^Ph. Several values of the magic number for 
Z have been predicted between 114 and 126, depending 
on the model used in the calculations, while that for N 
is 184 US H, IE 0- However, the major difficulty in 
the production of the SHE lies in the lack of knowledge 
about the reaction mechanism, especially about fusion 
mechanism of the massive systems. In the region of the 
heavy nuclei, the so-called fusion-hindrance is known to 
exist experimentallvP. lol. llo| . But no complete theoreti- 
cal description of the mechanism is yet given. As a con- 
sequence, it has been difficult to predict reliably which 
combination of target, projectile and bombarding energy 
is optimum to produce new elements. 

Let's assume the compound nucleus theory [lll|. though 
it might not be fully justified for such unstable systems. 
We can express the residue cross-section as 
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where J denotes a total spin of the system and Ec.m. 
the energy in the center of mass. E* is equal to Ec.m. + Q 
with the fusion Q- value. As usual P fusion ^^"^ P'surv ^-^^ 
the fusion and the survival probabilities for the spin J, 
respectively. The latter is obtained by the use of the sta- 
tistical disintegration code, which is newly dev elop ed by 
solving the time-dependent Bateman equation [l^. The 
former is the most unknown part. But recently, one of the 
present authors et al. have proposed the two-step model, 
where the fusion process is divided into two phases: the 



approaching phase and the formation phase. They cor- 
respond respectively to the passing over the Coulomb 
barrier and to the shape evolution toward the spherical 
compound nucleus, starting from the pear-shaped stick- 
ing configuration of the incident system. Thus, 



fusion sticking f orrnation 



(2) 



The calculation of the probabilities can be done by 
solving dynamics of the process in each step. The amal- 
gamated system is supposed to be exited, which mean 
the incident kinetic energy has been converted into the 
thermal energy. This conversion of kinetic energy would 
be done suddenly at the contact point or would start be- 
fore the top of the Coulomb barrier. We will assume the 
latter hypothesis, as it is well known that the dissipation 
observed in the deep inelastic collision (DIG) is under- 
stood by the viewpoint |l3j . With this assumption we 
calculate the sticking probability with a Langevin equa- 
tion including a frictional force and associated fluctuation 
with a time dependent temperature |l4l |. The formation 
probability is calculated by a multi-dimensional Langevin 
equation for the shape evolution with a constant temper- 
ature. The full description of this model has been given 
in the references [la, ll^ • 

We will first make a brief theoretical description of the 
two-step model. Then, we will check the ability of the 
model to reproduce known experimental data. Finally 
we will present predictions for the residue cross-sections 
for the reaction ™Zn -f^os Bi ^^78 _,_ 71q^ _,_208 

Pb ^278 ^_ 76g-g _,_208 ^283 ^ 

that have not yet been measured but should be in a near 
future. 



REMINDER OF THE TWO-STEP MODEL FOR 
FUSION 

We briefly recapitulate here the theoretical framework 
of the two-step model used to calculate the fusion prob- 
ability. 
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Approaching phase 

To describe the approaching phase, a classical descrip- 
tion of the relative motion between colliding nuclei is used 
with a frictional force and an associate fluctuation force. 
In the original treatment 13, the random force was 
not implemented, but is properly taken into account in 
the present model. 
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where fj, is the reduced mass and V is the sum of 
Coulomb and nuclear attractive potential. Ci{r) is the 
radial and the tangential friction, respectively. We will 
notice that here the rolling friction is neglected. Lq is 
the incident angular momentum and ^Lq so-called sliding 
limit. Ti{t) denotes a Gaussian random force with zero 
mean value. The equation0]is the dissipation-fluctuation 
theorem. In the case i — j = (j), factor is necessary 
for C(j,(r) in Eq. 0] In this modelisation the temperature 
increases as the energy is dissipated by the friction force. 

This modelisation is the one given by the surface fiction 
model (SFM)^3|i which has been used successfully to re- 
produce the main features of the deep inelastic collisions, 
expect the spin of the outgoing fragments which appear 
to require an inclusion of the rolling friction. Moreover 
the very different strengths between the radial and tan- 
gential frictional forces may not be physical, but may be 
a certain effective description. And of course, its applica- 
bility is limited to the energy region above the Coulomb 
barrier, due to its classical nature. Nevertheless, it is 
simple and has no free parameters, so we employ it for a 
description of the approaching phase for the realisation 
of the two-step model . 

One of the characteristic results of the model is the 
Gaussian distribution of the radial momentum at the 
touching point. The origin of this distribution comes 
from the Gaussian random force associated to the fric- 
tion force. The width of this distribution is consistent 
with the temperature determined by the internal energy 
transferred from the incident kinetic energy. This distri- 
bution is used as an initial condition for the subsequent 
dynamics of the shape evolution from a pear shape to- 
ward the spherical shape. 



Shape evolution 

The pear-shaped system formed by the contact of the 
nuclei involved in the collision is located outside of the 
conditional saddle point ^^H^]. As a consequence, we 
have to solve the shape evolution toward the spherical 
shape under a friction force associated with a random 



force. The dynamics of this process is again treated by 
the general multidimensional Langevin equation |2ll l22l 

Elm. 
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- -'-^-iirM~')^^w (5) 
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where V'^ denotes the liquid drop model (LDM) po- 
tential energy surface plus the rotational energy for spin 
J 0, 0|. The equation is again the dissipation- 
fluctuation theorem with temperature T'^ being assumed 
constant for the simplicity. The friction tensor is calcu- 
lated by using the so-called one-body model (OBM) , 
that is, a one-body wall-and-window formula with the 
two-center parametrisation of nuclear shape. 

The ratio between the number of trajectories reach- 
ing the spherical configuration and the total number of 
trajectories gives the formation probability. Mostly, the 
trajectories do not go over the ridge line but go back 
to re-separation. We calculate this probability for var- 
ious initial momenta at a given touching point. Then, 
the total formation probability is obtain by the convo- 
lution of those probabilities with the initial distribution 
of momenta obtained in the analyses of the approaching 
phase. 



COMPARISON WITH EXPERIMENT 

As the fusion probability is calculated as above, we 
will now calculate residue cross-sections according to Eq. 
^ In the scope of this paper we will only consider the 
\n reactions. The decay calculations are done for each 
partial cross-section for a given J by the code KEWPIE 
Il2j which has only one free intrinsic parameter that is 
the reduced friction coefficient (/3) used for the Kramers 
factor [13. This value is set to 13 = S.lQ-^is for aU the 
results presented in this publication in consistence with 
OBM. Here, the essential ambiguities are the masses and 
shell correction energies employed for SHE. To know their 
influence we have made the calculations employing three 
different tables of mass : MoUer & Nix, Koura, Myers & 
Swiatecki HIH's^. 

First, we wih take ^^Ni -f Pb ^^71 hq _^ system 
as an example for which more precise excitation function 
is recently measured at RIKEN, in addition to the pre- 
vious GSI data. Calculations are made with the fusion 
probabilities obtained above and the survival probabili- 
ties with KEWPIE code. The calculated peak positions 
in Ec.m. are not far from the data, though there are small 
difference among the mass tables. This is mainly due to 
the slight difference in the macroscopic part of the masses 
among the tables, which is consistent with the qualitative 
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FIG. 1: Residue cross-sections for systems that have been al- 
ready studied experimentally. The symbol represent experi- 
mental data with the error has in energy (thickness of the tar- 
get) and cross section (sensitivity of experimental device and 
cumulated dose). The line represent the results of calculation 
assuming different table of mass and shell correction energy. 



explanation by Swiatecki et al. n31j|. Absolute values of 
the peaks of In cross-sections are larger than the experi- 
ment by a few to several orders of magnitude depending 
on the mass tables. Nevertheless, it would be remarkable 
that the calculation of the present model result in rather 
good reproductions of the experiments without using any 
free parameter adjusted. For more precise reproduction 
and thus accurate predictions, however, we have to refine 
the model, in the calculation of PsUcking and/or Pform, or 
to find accurate values of the masses etc. But it would be 
practical to introduce a phenomenological reduction fac- 
tor, for example, a constant renormalization of the cross- 
sections. In view of a variety of predictions of the shell 
correction energy, we insteads introduce a phenomeno- 
logical factor for the shell correction energy which is to 
be used all through each mass table employed, keeping its 
own characteristic tendency predicted over the elements 
and the isotopes under investigation. An extreme sensi- 
tivity of the calculated peak height to the shell correction 
is another reason for the factor effective in the systematic 
analyses. 

In order to reproduce the peak heights of the experi- 
mental measurements for the reaction ^'^Ni-\-'^^^ Ph ^^71 
110 -I- In, the factor is 0.4 , 0.36 and 0.4 for the shell cor- 
rection energy of MoUer & Nix, Koura, Myers & Swiate- 



cki, respectively. 

In the Fig. ^ is shown the comparison between the 
calculations and the experiments |32L l33| . The energy of 
the experimental data plotted in this figure corresponds 
to the energy in the center of mass minus the energy loss 
at the middle of the target. The associated error bars in 
energy correspond to the energy losses in the target. 

In the upper left panels of the Fig. ^we notice that 
the position of the peak in Ec.m. is well reproduced, es- 
pecially for the tables of MoUer & Nix and Myers & 
Swiatecki. For a justification of the modification fac- 
tor for of the shell correction we have made calcula- 
tions for the systems ^^Ni -^^09 bi ^272 m and 
™Zn -1-208 P5 ^277 + In, keeping the same value of 
the multiplication factor for the shell correction energy. 
The lower panels of the figure Q show an overall (cal- 
culated curve within the experimental error bar) good 
agreement between the experimental measurements and 
the theory. Moreover, the global phenomenological de- 
creasing trend of the peak heights of the cross-sections is 
well reproduced, assuming the MoUer & Nix and Myers 
& Swiatecki. For those two mass tables we can even say 
that a remarkably good agreement between the experi- 
ment and theory is obtained. We should remind that the 
same model has been already applied to the hot fusion 
path {^^Ca -\- actinide reaction) and resulted in a fairly 
good reproduction of the data 0, 0| . 



PREDICTIONS OF THE RESIDUE CROSS 
SECTION FOR ELEMENTS Z=113 AND Z=114 

As we are able to reproduce very well the residue cross- 
section for the above known systems we can extend this 
recipe to the systems where measurements are not yet 
made. We will focus on the yet undiscovered element 
Z=:113 and a new isotope of the Z = 114 produced by 
the cold fusion path. In order to calculate residue cross- 
sections we will just extend the calculations with the 
same factor for the shell correction energy as that we 
have set for the ^^Ni -l-^os Ph ^^71 no + In reaction. 

For the production of the Z=113 the reaction ™Zn-|-20^ 
Bi ^278 _^ ^i^g reaction '^^Ga ^"^^^ Pb ^278 

113 + In is under consideration. As we see on the upper 
panels of Fig. |21 both reactions are predicted to lead to 
roughly the same residue cross-sections for all the mass 
predictions employed. This similarity come from the fact 
that on one hand the fusion barrier is unfavourable for 
the reaction '^^Ga+'^^^Ph ^27i 113 + 1^ but on the other 
hand the high Q-value of this reaction leads to lower ex- 
citation energy that is favourable for the emission of neu- 
tron in the competition between fission and evaporation. 

Moreover, the results are in agreement with the global 
phenomenological trend of a decreasing of the In residue 
cross-sections by roughly a factor 3 for each increasing of 
one charge unit in the SHE. 
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FIG. 2: In Residue cross-section for system that have not been 
yet measured. The line represent the results of calculation 
assuming different table of mass and shell correction energy. 



The production of the element Z—llA is the other case 
that we will predict. This element has been already pro- 
duced by the hot fusion path |34l| but the experiment 
with the reaction "^^Ge+^^^Pb 114+ In will be the 

first measurement of the production of Z = 114 by the 
way of a In reaction (cold fusion path). The results of 
the calculation are shown in the lower panel of Fig. |21 
The predicted value is around a few tenth of pico-barn. 
To measure the residue with a so low cross-section is a 
challenge in experiment. 



CONCLUSION 

As we are able to reproduce rather well both the po- 
sitions and the absolute values of the peaks of measured 
In residue cross-sections by introducing only one free 
parameter, one reduction factor for the predicted shell 
correction energies, the two step model appears to be 
quite promising. Based on the success we have made 
the prediction on the In cross-sections in ™Zn -|-209 
Bi ^278 ;li3 + In, ^^Ga +^°^ Pb ^^78 ^^13 + i„ and 
™Ge -1-2°^ Pb ^283 ^ 2n without adjusting anything 
more. This is the first predictions of the excitation func- 
tion of residue cross section for Z=113 and 114 by the 
calculation with the dynamical model of reactions. Of 



course, the ultimate test of the model is the future ex- 
periment performed on the above systems. The model is 
expected to provide a reliable guideline for future exper- 
iments for productions of super-heavy elements. Hope- 
fully the present model permits have an overview of the 
most suitable experiments and win a lot of experimental 
time. 
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